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Abstract 

We discuss the renormalization of a BRST and anti-BRST invariant com- 
posite operator of mass dimension 2 in Yang-Mills theory with the general 
BRST and anti-BRST invariant gauge fixing term of the Lorentz type. The 
interest of this study stems from a recent claim that the non-vanishing vacuum 
condensate of the composite operator in question can be an origin of mass gap 
and quark confinement in any manifestly covariant gauge, as proposed by one of 
the authors. First, we obtain the renormalization group flow of the Yang-Mills 
theory. Next, we show the multiplicative renormalizability of the composite 
operator and that the BRST and anti-BRST invariance of the bare composite 
operator is preserved under the renormalization. Third, we perform the op- 
erator product expansion of the gluon and ghost propagators and obtain the 
Wilson coefficient corresponding to the vacuum condensate of mass dimension 
2. Finally, we discuss the connection of this work with the previous works and 
argue the physical implications of the obtained results. 
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1 Introduction 



It is still a challenging and unsolved problem to prove quark confinement in the 
framework of quantum chromo dynamics (QCD). A very beginning question in de- 
riving quark confinement is in what sense quark is confined? A simple criterion of 
quark confinement which has been widely used so far is the area law decay of the 
Wilson loop (defined by the vacuum expectation value of the Wilson loop operator). 
The area law implies the presence of a linear piece ar proportional to the interquark 
distance r in the static interquark potential V{r). The dual superconductivity of 
QCD vacuum [jl| is one of the most promising mechanisms of quark confinement in 
compatible with this picture. However, it is well known that this criterion is not so 
useful in the presence of dynamical matter, since the interquark force is screened by a 
pair of quark and anti-quark created from the vacuum and the linear piece no longer 
appears in the potential. 

In the previous paper 0, one of the authors (K.-I. K.) has proposed a non- 
vanishing vacuum condensate {O) of mass dimension 2 as the origin of mass gap and 
quark confinement in Yang-Mills theory. The proposed composite operator of mass 



dimension 2 is given by 

0:= J^/rf-xtr 



-Af,{x) ■ A^{x) + aiQ{x) ■ Q{x) 



where A^j, is the gauge field, C (C) is the ghost (anti-ghost) field and fi^^^ denotes the 
volume of the D-dimensional spacetime. It has been shown that the composite 
operator O is invariant under the Becchi-Rouet-Stora-Tyutin (BRST) and anti- 
BRST 0] transformations in the manifestly Lorentz covariant gauge, especially in 
the most general^ Lorentz gauge [|5|-p!0[1 and the Maximal Abelian (MA) gauge [pH] - 



In ( |1.1|) , the trace is taken over the broken generators of the Lie algebra Q of 
the original group G when the original gauge group G is broken to by a local 
gauge fixing condition chosen, i.e, G itself in the Lorentz gauge and G/H in the 
MA gauge corresponding to the maximal torus group H of G. Especially, in the 
limit a — > (which we call the Landau gauge), the composite operator reduces to 
O = (Q^^^)"^ J d^x tr |yi^(x) ■ A^{x) and hence it becomes gauge invariant, since 
the contribution from the ghost and anti-ghost disappears. The vacuum condensate 
includes the ghost condensation proposed in the MA gauge and reduces to the 



gluon condensation proposed recently by several authors [piH 24|, see also p5i|26[ . 

The physical implication of the existence of such a condensate {O) has been argued 
based on the operator product expansion (OPE) of the gluon and ghost propagators 
(2-point functions) and the vertex function (3-point function) |^^,^. However, 
the actual calculation has been performed within the tree level. 

In order for such a proposal to be meaningful, it is very indispensable to show that 
the whole strategy to derive quark confinement based on the novel vacuum condensate 



"'^The precise definition of 'the most general' is stated later in the text. Roughly speaking, the 
most general Lorentz gauge is obtained by imposing both the BRST and anti-BRST invariance for 
the gauge fixing term which corresponds to the Lorentz gauge d^A^{x) — 0. The resulting gauge 
fixing term has two parameters. The conventional Lorentz gauge is obtained as a special choice of 
the parameters. 
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survives the renormalization. In view of this, we focus on the renormahzation of the 



composite operator (|1 . 1| ) . The main purpose of this paper is to examine whether or 
not the composite operator in the integrand of O is renormahzable. In addition, we 
must clarify the meaning of the BRST and anti-BRST symmetry in the renormalized 
theory. We examine whether or not the renormalized composite operator is 
invariant under the renormalized BRST and anti-BRST transformation. If this is the 
case, the proposed composite operator of mass dimension 2 and the corresponding 
vacuum condensate can have a definite physical meaning. The analysis of this paper 
is restricted to the most general Lorentz gauge fixing, since the analysis of the MA 
gauge is more involved and hence the result is to be reported in a separate paper p7| . 

In the most general Lorentz gauge, the multiplicative renormalizability of the 
Yang-Mills theory has been worked out by Baulieu and Thierry-Mieg |^ by making use 
of the Slavnov- Taylor identities characterizing the BRST and anti-BRST invariance of 
the theory (See textbooks and reviews, e.g., pB|-P^ ). In the course of renormalizing 
the composite operator, however, there is a subtle problem of the operator mixing. 
In order to discuss the renormalization of a composite operator, we must take into 
account all the contributions coming from all the other composite operators of the 
same mass dimension and the same symmetry property. In the OPE, the Wilson 
coefficient corresponding to arbitrary vacuum condensate can be calculated by the 
perturbation theory. In the usual Lorentz gauge, the Wilson coefficient associated 
with the ghost condensate (C ■ C) in the OPE of the propagator vanishes identically 
due to a special property of the 3-point gluon-ghost-anti-ghost vertex as pointed out 
in In the most general Lorentz gauge 1^,0, however, we show in this paper 

that the operator mixing between two composite operators, ^A^ ■ A^^ and iG ■ C, of 
mass dimension 2 does exist in general due to the presence of four-ghost interaction 
(except for the case which is reduced to the conventional Lorentz gauge). We explicitly 
calculate the matrix of renormalization factors of the composite operator in the one- 
loop level. 

For the Landau gauge, the vacuum condensate of mass dimension 2 in Yang- 
Mills theory is nothing but the gluon pair condensation. A possibility of gluon pair 
condensation was already suggested from the existence of the tachyon pole in the 
two gluon channel by solving approximately the Bethe-Salpeter equation, see e.g. |^| 



and [0. A gluon pair can be identified as a Cooper pair, that is a bound state 
caused by the attractive force. Hence the gluon condensation is regarded as the Bose 
condensation of the gluon with spin 1. A remarkable point of our treatment different 
from the previous one is to retain the manifest Lorentz covariance and gauge (or BRST 
and anti-BRST) invariance. Hence the introduction of ghost field is indispensable in 
this approach. It is important to clarify how the inclusion of the ghost infiuences the 
dynamics of gluon to recover the gauge invariance. This paper is a preliminary work 
toward the complete understanding of this problem. 

Another purpose of this paper is to point out that the composite operator discussed 
above has the analog in the Abelian gauge theory, especially, quantum electrodynam- 
ics (QED). This suggests that a confinement phase can exist even in QED probably 
in the strong coupling region [39-42]. In QED, the vacuum condensate in question is 
reduced in the Landau gauge to the photon pairing. The photon pairing has also been 
suggested long ago by solving the Cooper equation, see |H3lH|. From quite a different 
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viewpoint, one of the authors |^ discussed the existence of a confinement phase in 
QED based on the total QED Lagrangian with the BRST and anti-BRST invariant 
gauge fixing term which is identical to the usual Lagrangian in the Lorentz gauge up 
to a total derivative term. An advantage of rewriting the the gauge fixing part of the 
Lagrangian into the BRST and anti-BRST exact form is that the hidden supersym- 
metry becomes manifest and that the gauge-fixing part in four spacetime dimensions 
is reduced to the 0(2) non-linear sigma model in two spacetime dimensions owing to 
the Parisi- Sour las dimensional reduction.^ In view of this, the ghost is indispensable 
in this approach even for the Abelian gauge theory where the ghost decouples and is 
usually considered to be unnecessary. In the analysis of quark confinement, it is most 
important to understand the origin of the scale or the mechanism of mass generation 
which was not so clear in the previous treatments. The detailed analysis of this issue 
will be reported in a subsequent paper. 

This paper is organized as follows. In section 2, we summarize the BRST and 
anti-BRST transformations and their properties which are necessary in the following 
analyses. 

In section 3, we examine how the renormalization in QED is performed so as to 
preserve the BRST and anti-BRST symmetry. This section is a preliminary step for 
dealing with the non- Abelian gauge theory in the subsequent sections. 

In section 4, we consider the most general Lagrangian of the Yang-Mills theory 
which has manifest Lorentz covariance, global gauge invariance and BRST and anti- 
BRST symmetry. The gauge fixing term contains two gauge fixing parameters. We 
give Feynman rules of this theory and calculate the renormalization constants in the 
one-loop level. Although some materials in this section are a reconfirmation of the 
results obtained by Baulieu and Thierry-Mieg [§, it is necessary to make this paper 
self-contained and to give basic ingredients in the subsequent sections. 

In section 5, we obtain the renormalization group fiow in the parameter space of 
the theory. To one-loop order, e specify the location of the fixed points and obtain 
the equation of the lines of connecting the fixed points. 

In section 6, we discuss the main subject of this paper; the renormalization of 
the composite operator O of mass dimension 2. First, we show when the composite 
operator O is both BRST and anti-BRST invariant. Next, we evaluate the renormal- 
ization of O by taking into account the mixing of the operators with the same mass 
dimension and the same symmetry. To the best of our knowledge, the renormalization 
of the composite operator of mass dimension 2 has not been fully discussed except for 



a special case, i.e., the Landau gauge in the conventional Lorentz gauge fixing [p^ . 

In section 7, we perform the operator product expansion of the gluon and ghost 
propagators and obtain the Wilson coefficient associated with the vacuum condensates 
in question. 

In the final section, we give the conclusion of this paper and discuss the future 
directions of our research. In Appendix, we give some of the calculations omitted in 
the text. 



^ This formulation has been appUed to QED at finite temperature and a new confining phase is 
claimed to exist, see Mq] and references therein. 
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2 BRST and anti-BRST transformation 

We consider the general non-Abelian gauge theory with a gauge group G. In the 
following we use the notation: 

F-G:=F^G^, F^:=F-F, {F x G)^ := f^^^ F^G^ , (2.1) 

where f^^'^ are the structure constants of the Lie algebra S of the gauge group G. 
For the non-Abelian gauge theory, we define the BRST transformation by 

dsA^ix) = D^[A]eix) := d^eix) + g{A^{x) x e{x)), (2.2a) 
SbG{x) = -^g{e{x) X e(x)), (2.2b) 

SBGix) = i'Bix), (2.2c) 
dB'B{x) = 0, (2.2d) 

where yi^, "B, C and C are the non-Abelian gauge field, the Nakanishi-Lautrup (NL) 
auxiliary field, the Faddeev-Popov (FP) ghost and anti-ghost fields respectively. An- 
other BRST transformation, i.e., anti-BRST transformation is defined by 

SBA^ix) = D^,[A]eix) := d^eix) + g{A^{x) x e(a;)), (2.3a) 
dBGix) = -i(?(e(x) X e(a;)), (2.3b) 

SBe{x) = iS(x), (2.3c) 
Sb'B{x) = 0, (2.3d) 

whereQ S is defined by 

'B{x) = -'B{x) + ig{e{x) x e(x)). (2.5) 

The BRST and anti-BRST transformations are nilpotent and they anti-commute: 

SbSb = 0, 6BdB = 0, SbSb + SbSb = 0. (2.6) 

For the Abelian gauge theory, the BRST transformation reads 

SBa,iix) = <9^C(x), (2.7a) 

5bC{x) = 0, (2.7b) 

6bC{x) = iB{x), (2.7c) 

SbB{x) = 0, (2.7d) 



^ The last transformation is equivalent to 

Sb'B{x) ^ -gGix) X -Bix). (2.4) 
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where A^, B, C and C are the Abehan gauge field, the NL auxihary field, the FP 
ghost and anti-ghost fields respectively. The anti-BRST transformation is reduced to 



5Ba^,{x) = d^C{x), (2.8a) 

6bC{x) = 0, (2.8b) 

SbC{x) = iB{x), (2.8c) 

dBB{x) = 0, (2.8d) 



where B is defined by 



B{x) = -B{x). (2.9) 



3 QED in the Lorentz gauge 

As a warming-up problem, we consider the quantum electrodynamics (QED). As is 
well known, the total Lagrangian of QED is given by 
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^QED = -^rf^u + i^{il^d, - m)^ - e^i^i,a^ + £gf+fp, (3.1) 

with a gauge-fixing (GF) plus FP ghost term £jgf+fp- The explicit form of the 
GF-I-FP term depends on the gauge chosen. In this paper we adopt the most familiar 
covariant gauge, i.e., the Lorentz gauge 

d^a^ = 0. (3.2) 

Therefore, the GF+FP term is given by 

'Cgf+fp = -^^B (Cd'^a^ + ^CB^ = Bd^a^ + ^B^ + iCd^d^C. (3.3) 

Although the ghost and anti-ghost fields are free and decouple from other fields, we 
have included them to study the relationship with the non-Abelian case which will 
be discussed in the next section. 

As pointed out in |3^, the GF-I-FP term (|3^ ) is rewritten into the BRST and 
anti-BRST exact form. 



Lgf+fp = iSBSB ( ^a^a^ + ^iCC ) . (3.4) 



In fact, this is cast into the form, 



^GF+FP = i^B (^{dBa^)af, - '^iCSBC^ 



i5B{d^Ca^-^CB) , (3.5) 



which agrees with (|3.3| ) up to a total-derivative term. 
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If the NL field B is eliminated by performing the functional integration or by 
making use of the equation of motion, then we obtain 

'C^GF+FP = -^{d'^a^y + iCd'^d^C. (3.6) 
za 

The on-shell BRST transformation is given by 

SEdfiix) = df,C{x), (3.7a) 

6bC{x) = 0, (3.7b) 

SBCix) = --d^a^ix), (3.7c) 
a 

while the on-shell anti-BRST transformation is 

^Ba,j.{x) ^ du_C{x), (3.8a) 

6bC(x) = 0, (3.8b) 

6bC(x) = +-d^aJx). (3.8c) 
a 

The GF+FP Lagrangian ^qf+fp ^^<i total Lagrangian ^q^y) "^i^^ '^gf+fp ^^P" 
arately invariant under the on-shell BRST and on-shell anti-BRST transformations. 
The nilpotency of the on-shell BRST and anti-BRST transformation is realized only 
when the equation of motion for the ghost and anti-ghost field is used, since 

(SBfa^ix) = 0, (3.9a) 

(5b)'C(x) = 0, (3.9b) 

{SByC{x) = -'-d^'d^Cix), (3.9c) 

and 

\2 



{5BYa^{x) = 0, (3.10a) 



{5bYC{x) = +-d^'d^C{x), (3.10b) 
a 

{5BfC{x) = 0. (3.10c) 
Moreover, we obtain the similar result for the anti-commutability: 

{6b6b + ^b«5b)s(x) = 0, (3.11a) 

{8b5b + 5bSb)C{x) = --d^d^C{x), (3.11b) 

a 

{SbSb + dBdB)C{x) = +-d''di,C{x). (3.11c) 

a 

Now we define the composite operator O of mass dimension 2 by 

C := ^ / d^'xQix), Q{x) := ^a^(x)a^(x) + aiC{x)C{x). (3.12) 



6 



This composite operator is BRST and anti-BRST invariant, since 

6bQ{x) = d''{a^{x)C{x)), dBQ{x) = d>'{a^{x)C{x)). (3.13) 

We consider the renormahzation of the composite operator Q. The Abehan case is 
very simple due to the trivial renormahzation factors , Zee for the composite field 
^a^a^ and iCC. Therefore, we have only to take into account the renormahzation 
factor of the fundamental field, a^, C, C and the gauge fixing parameter a. QED 
is known to be multiplicatively renormalizable in the sense that the divergences are 
absorbed by introducing the renormahzation factors in the following way. 

V. = Z2^/V, (3.14) 

a, = Z\!\l, (3.15) 

C = ZcC^, C = ZcC^, (3.16) 

{B = Z^^^^B^), (3.17) 

m = Z^Z^^m^, (3.18) 

a^Zo^a^, (3.19) 

e = ZiZ-^Z^^^^e^. (3.20) 

The renormahzation factors are not independent to each other. In fact, the coupling 
constant is renormalized as 

e = Zg-'/'e^, (3.21) 
as a consequence of the Ward relation: 

Zi = Z2. (3.22) 
Moreover, the Ward-Takahashi identity yields 

Za = Zs. (3.23) 

The result of perturbative renormahzation in QED is well known and can be seen in 
the text books. A result: 



(3.24) 



means that both ghost and anti-ghost arc free and receive no renormahzation in the 
perturbation theory (This is not the case in the non-Abelian case). Consequently, we 
arrive at the result that the composite operator is renormalized as 

Q = ZaQ^, := \a^{x)a^\x) + o^iC\x)C\x). (3.25) 

Therefore, the BRST invariant combination of two composite operators with mass 
dimension 2 is preserved under the renormahzation. 

In view of the above results, the renormahzed BRST transformation is defined by 

dl^zy^B, 61 = Z]IHb. (3.26) 
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This is shown as follows. The Noether current of the BRST symmetry is obtained as 

= Bd^'C - d^'BC - duiF^C). (3.27) 

The Noether charge, i.e. the BRST charge as the generator of the BRST trans- 
formation 



is given by 



[zAQs,$(x)] = A5b$(x), (3.28) 



Qb = j d^xJl = J Sx[Bd^C - d^BC]. (3.29) 



In the similar way, the anti-BRST charge Qb can also be defined as the Noether 
charge for the anti-BRST transformation. Therefore we can define the renormalized 
BRST charge Q% as 

Ql = zI^^Qb = I - a°5^C^]. (3.30) 



This ensures the renormalization of the BRST transformation ( ^.26] ). The renormal- 



ized BRST transformation for the renormalized field has the same form as the bare 
BRST transformation for the bare field. Thus, the composite operator Q is a BRST 
invariant and multiplicatively renormalizable operator for arbitrary gauge parameter 
a. The renormalized GF+FP term has the same form as the bare one: 

i^GF+FP = iS^di f ^ + ^zC'cA . (3.31) 



4 Yang-Mills theory in the most general Lorentz 
gauge 

4.1 Lagrangian 

We consider the most general quantum Lagrangian density that is a local function 
of the fields, A^, "B^, Q^, C"^ and satisfies the conditions: The Lagrangian is (1) 
of mass dimension 4, (2) Lorentz invariant, (3a) BRST invariant, (3b) anti-BRST 
invariant, (4) Hermitian, (5) of zero ghost number, (6) global gauge invariant, and 
the theory with this Lagrangian is (7) (multiplicative) renormalizable. Here it is 
implicitly assumed that the Lagrangian is written as the polynomial of the fields, 
and that there are no higher derivative terms, since there is no intrinsic mass scale 
in the Yang-Mills theory. It should be remarked that we have imposed BRST and 
anti-BRST invariance instead of gauge invariance (we do not require gauge invariance 
for the Lagrangian). Such a Lagrangian was given by Baulieu and Thierry- Mieg 
as 



+ i5^5j, (a^A^ ■ A^ + ■ e) + -S ■ S, (4.1) 
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where ai{i = 1, 2, 3, 4) is an arbitrary constant, and and Sb are the BRST and anti- 
BRST transformations. The first term is the Yang-Mills Lagrangian, the second term 
is the topological term which is not discussed in this paper and omitted hereafter. 
The first and the second terms are gauge invariant. On the other hand, the third 
and the fourth terms are identified with the GF and FP term, since they break the 
gauge invariance of the Lagrangian. After the rescaling of the parameters and the 
field redefinitions, we can cast the total Lagrangian of the Yang-Mills theory into the 
form, 

^YM = ~^^fJ.!^ ' ^^'^ + ^GF+FP, (4.2) 

with the GF+FP term [|-[10 



£gf+fp = iSbSb (^^A^ . yiM - ■ e) + ■ :B (4.3) 

= -i8b [-d,Q .A^ + -Q-'B- -age ■ (C x 6) j +^2-3. (4.4) 

The final term is allowed for the renormalizability of the total Lagrangian and is 
written in either BRST exact or anti-BRST exact form, 

S-S = -i(5B(e-S)=i^B(e-3). (4.5) 



However, the GF+FP term ( ^.4] ) is simultaneously BRST and anti-BRST exact, i.e., 
5b5b(*), only if a' = 0. If we impose one more condition, e.g., the FP ghost conju- 
gation invariance, 

e^^±e^, e^^^e^, s^-^-s^, ^a^^a^), (4.6) 

the second term of ( [4. 4] ) is excluded, namely, only the choice a' = is allowed. 
By performing the BRST and anti-BRST transformations, we obtain 

^GF+FP = . s - |*^?(e X e) ■ B + s ■ d,A^ 

+ ie ■ d^D^'[A]e + -g\e x e) ■ (e x e) (4.7) 

8 

= ^ ^ _ i^^^g X e) ■ B + S ■ d^A^ 

+ ie ■ d^'D''[A]e + ^g\ie x e) ■ {ie x e). (4.8) 

The GF+FP term includes the ghost self-interaction where the strength is propor- 
tional to the parameter a. 

When a = 0, this theory reduces to the usual Yang-Mills theory in the Lorentz 
type gauge fixing with the gauge fixing parameter a': 

a' 

£gf+fp = ■ S + 3 ■ d^A^" + ie ■ d^'D''[A]e. (4.9) 
This is consistent with the FP prescription. 
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When a 7^ 0, there exists a quartic ghost interaction which can not be implemented 
by the usual FP prescription. Therefore we must go beyond the FP prescription. The 
GF+FP term is further rewritten as 

^GF+FP = - ^(d^^.f + (1 - O^e • d,'D^'[A]e + e^e • 'D^[A]d,e 
+ ^Ae(i-0/(^exe)-(zexe) 

+ ^ (S + X-'d^'A, - ^zg{e X e))' (4.10) 

= - ^{d^'A^f + te ■ d,,d^Q - (1 - OgiA^' ■ {d^Q x C) 

+ igiA^ ■ (e X d^Q) + iAe(l - i)g\i<^ x C) ■ {iQ x C) 

+ ^ (S + A^^S^'yi^ - iig{Q X e))' , (4.11) 
where we have defined the two parameter^ 

X:=a + a\ ^ := ^ = ^- (4-12) 

In this form, it is easy to eliminate the Nakanishi-Lautrup field B. We call the gauge 
( |4.11 ) the most general Lorentz gauge hereafter. 



4.2 Feynman rules 

We obtain the following Feynman rules for the Yang-Mills theory of the Lagrangian 
( ^12] ) with ( |4.11| ) where the NL field is eliminated. 



4.2.1 Propagators 

(a) gluon propagator: 



y4,/ixA/WWW»-B,z/_ Y^AB ^ 



g^y - (1 - A) — 
pi 



^AB 



(4.13) 



(b) ghost propagator: 



A ^— B = ,qAb ^ _ 1 ^AB^ ^4_^4) 



4.2.2 Three-point vertices 

(c) Three-gluon vertex: 



A. 

^ ^/^^^ [(g - '^)/.^P<x + (r - p)p^^^ + (p - g).^^p] . (4.15) 

p %!5\r 



^ The parameters a, a',A, ^ in this paper corresponds respectively to the Ac,Afc,A, a in |^ and 
a, a' , A, a/2 in m. 
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(d) Gluon-ghost-anti-ghost vertex: 



•^l^rJiN^^ =igf^''''[ap-q)-pf. (4.16) 



B\ 



4.2.3 Four-point vertices 

(e) Four-gluon vertex: 



ECDt 




(4.17) 

where I^u,pcT ■= {gtipQua - gi,agup)/'2 • 
(f) Four-ghost vertex: 

X = - 0/ (r^^r^^ - f^^^fBCE^ . (4.18) 

/ \ 

4.3 Multiplicative renormalization 

It has been proved by Bauheu and Thierry-Mieg based on the mathematical induc- 
tion that the Yang-Mills theory in the most general Lorentz gauge ( [4.11| ) is multiplica- 
tively renormalizable. We introduce the renormalization constant (or renormalization 
factor) for the field: 

= zK^Al, e = Z]i^Q^, e = Z]i^t^, S = Z^/^S^ = ZcZ^^'^'B^, (4.19) 
and for the parameters: 

A = ZaAr, i = Z^U. 9 = Zggn. (4.20) 
By substituting (^.19]) and (^.20|) into the bare Lagrangian, we obtain the total 



Lagrangian written in terms of the renormalized fields, renormalized parameters and 
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the renormalization factors: 



f tot 



1 



Rn2 



1 



- (1 - Z^U)Z^zi^ZcgK^A^'' ■ (9^e^ X e^) 

+ z^z.z^ZcUgi.iA^'' ■ (e^ x a^e^) 

+ ]^z^z^zlzl\^i^{\ - z^i^)gl{i<z'' X e^) • (ze^ x e^) 

+ ^Z, (ZcZ^'/'S^ + Zl^Z^X^d^Al - Z^Z.ZcUig^^'' x 



(4.21) 



The total Lagrangian ( 4.21|) is decomposed into a renormahzation-factor independent 
part ^YM^ ^^"i the remaining part ^Cym^ 



r tot 
'^YM 



tot R 
YM 



p tot c 
■*-'YM 



r tot R _|_ P tot c 
■ '^YM ~r 



-■YM ; 



(4.22a) 



--{dX-duAl^g^Al^A'^f 



2A 



R 



- (1 - U)g^iA^'^ ■ (9^e^ X e^) + Ug^iA^"" ■ (C^ x 9^6 

+ - iK)gl{}^'' X e^) ■ (^e^ x e^) 



R^ 



(|42l| ) - (i]22D. 



(4.22b) 
(4.22c) 



Here ^.'ym^ is obtained by setting all renormalization factors Z = 1 in ( [4.21|) and 
hence it is written in terms of the renormalized fields and renormalized parameters 
and has the same form as the bare Lagrangian -C-ym) while ^^Ia is the counterterm 
defined by the difference £ym ~ '^ym^- 



4.3.1 Renormalization of two-point functions 

First, we calculate the vacuum polarization function of the gluon. To the order g^^ 
there are three Feynman diagrams, see (al), (a2) and (a3) in Fig.|l|. 



fall 



vAA/( 



(a2) jJv/Wu 




Figure 1: Vacuum polarization of the gluon. 
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As a gauge-invariant regularization, we adopt the dimensional regularization. 
Then we obtain the following result (e := 2 — Z^/2). 



(al) = C2{G)5 



(a2) = -C2{G)5 



(4vr) 
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e(i-0 



(a3) 



2 
0, 



(4 



TT 
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T 1 



(4.23a) 

yg^g^ + (1 - \)iq^g^,u - q^,qu)^ , (4.23b) 

(4.23c) 



where C2 = C2{G) is the quadratic Casimir operator in the adjoint representation of 
the gauge group G defined by S^^C2{G) = f^CD ^bcd ^ Hence the counterterms 
and 5l for the transverse and longitudinal part of the vacuum polarization tensor are 
determined so as to satisfy the relation, 



(al) + (a2) + (a3) - i5T:{q^g,, - g^g.)(5^^ - ^yg^g.^^^ 



0, 



which yields the result: 

A\ {g^i-^fG2{G) 
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(47r)^ 



-Ae(l-0 



{g^i'^fC2{G) 
(47r)2 e ' 



On the other hand, the relationship 



Za - 1 



'(1) 



+ 



7(1) 



+ 



(4.24) 



(4.25) 



(4.26) 



must hold for the multiplicative renormalizability where we have defined the renor- 
malization factor Z order by order of the loop expansion, Z = 1 + Z'^^'^ + Z^'^'^ + • • • . 
Thus we obtain the renormalization factors: 



13 X\{gii~^fC2{G) 



and 



St — 5t 



13 
~6 



{gfi-rC2{G) 
(47r)2 e ■ 



(4.27) 



(4.28) 



Note that 6t and hence Za is the same as in the FP case where the four ghost 
interaction does not exist. When ^ 7^ 0, 1, however, we find that 5l 7^ or equiva- 
lently Z^ 7^ Zx- On the contrary to the FP case, the longitudinal part of the gluon 
propagator must be renormalized in this case. 

Next, the vacuum polarization function of the ghost is calculated in the similar 
way. To the order g"^, there are two Feynman diagrams, see (bl) and (b2) in Fig.^. 
The explicit calculation shows that 



(bl) 
(b2) 



igf^'^YC2iG) 



4 



{Any 



0. 



(4.29a) 
(4.29b) 
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Figure 2: Vacuum polarization of the ghost 



The counterterm 6c is determined from 



(bl) + (b2)-p^5^^5c = 0. 



(4.30) 



Hence the counterterm Sc = Zq — 1 = Z^^ + ■ ■ ■ is equal to the renormalization 



constant Z(^^: 



4^ = Sc 



4 {Any 



(4.31) 



This is again the same as in the FP case. 

4.3.2 Renormalization of three-point function 

We consider the renormalization of three-point vertex. For example, the Feynman 
diagrams for the radiative correction of the gluon-ghost-anti-ghost vertex to one-loop 
order is given in Fig.|[ 



(cl) 




(c2) ^ 



(c3) / 



Figure 3: Radiative corrections for the gluon-ghost-anti-ghost vertex. 



If we write the counterterm for the gluon-ghost-anti-ghost vertex function as 



we find the renormalization factors are related as 



^Acc - ZcZA^g^i 



\z,z, - 1 = z^.^) + iz« + Z« + Zf + 



b\cc = ZcZiz, - 1 = Z« + ^4^) + + ■ ■ ■ . 
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(4.32) 



(4.33) 
(4.34) 
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At p = q, the respective diagram is calculated as 



(c3),=, = 0. 

By substituting (^ISSa]) , (|]35|) and (|]35|) into 

(cl)^^^ + (c2)^^^ + (c3)^^^ - tgf^''^5lccP, = 0, 

it follows that 

Hence the renormalization factor is obtained as 

<? Acc c 2 ^ 6 (47r)2 e 

At p = 0, the respective diagram is calculated as 



. i 13 



By substituting (W^) , (HH) and (|]39|) into 

(cl)^^o + (c2)p^o + (c3)^^o - tgf^''^U5],ccQ, = 0, 

it follows that 

{gfi-^fC2{G) 
(47r)2 e ■ 



-A(l-OI?-i)-iA; 



Then we obtain 



'(1) _ Xl _ _ X2 _ 

'5 ~ "ACC "acc 



2) (47r)2 e 

Accordingly, the renormalization constants of a and a' are obtained as 

7(1) =(]1_'^\ i9f^~'?C2iG) 
6 a) (47r)2 e ' 



and 



13 a + a'\ {giJ~'fC2{G) 
6" 2~j (47r)2 ~' 
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5 Renormalization group flow and flxed points 



Using the above result, the renormahzation group (RG) functions are obtained as 
follows. The /3-function is obtained as 



(5.1) 



It turns out that the /^-function does not depend on the gauge parameters A and ^: 



dgR 



1 11 



(9/i IQit"^ 3 
Similarly, we obtain the RG functions: 

d_ 

and 



C2{G)gl 



7S := f^^^R = ^XkUCk - 1) (^R - 



d 

7a := At^-^R = 2Ar 



^ -^R^r(1 — ^R 

o 2 



l) (47r)2^^ 
C2(G) 



(47r) 



-9k- 



(5.2) 



(5.3) 



(5.4) 



The RG flow in three-dimensional parameter space (^, A, g) is determined by solv- 
ing simultaneous differential equations: 



^dil - ^^^^^ " " 2 j l4^' 



dX 

/i— = 2A 



13 A 
Y ~ 2 



(9/i 

dg _ 11 C2{G)g' 
9/i 



+ Ae(i-o 



g2(g)^^ 

(47r)2 



3 (47r) 



(5.5a) 
(5.5b) 
(5.5c) 



where we have omitted the subscript R for the renormalized quantity. 
As is well known, the equation ( |5.5c| ) is solved exactly. 



1 



1 + 



22 C2(G) ^2 
3 (47r)2 



g^if^o) In 



_M 22 CijG) 
3 (47r)2 



MO 



In 



(5.6) 



where we have used the boundary condition g{fio) = oo at /io = Aqcd- The remaining 
two equations (|5.5a|) and ( |5.5b|) can not be solved exactly. 



5.1 Fixed points 

First, we obtain the fixed point of the RG. Note that the derivative ^A*^ in ( 5.5a| ), 
( p.5b| ) is rewritten as 



1 d 22C2(G), /i d 22C2(G) d 



^2^a/x 3 (47r)2 Aqcd dl^ 3 (47r)2 51nln- 



(5.7) 



^QCD 
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Then the fixed point (to one-loop order) is obtained by solving the algebraic equation 
simultaneously: 



Ae(e-l)(e-^) =0, A 



0. (5. 



We find one line of fixed points and three isolated fixed points in the A) plane or 
equivalently four isolated fixed points in the (a, a') plane: 

A. The line of fixed points: A = 0,^ G R corresponds to an isolated fixed point 
{a, a') = (0,0). 

B. (^, A) = (i, f ) corresponds to (a, a') = (f , 0). 

C. (^, A) = (0, f ) corresponds to {a, a') = (0, f ). 

D. (^, A) = (1, f ) corresponds to (a, a') = (f , -f ). 

If the two parameters ^, A are set equal to one of the fixed points, the theory remains 
forever on the fixed. If the system starts from other points and the scale /i is decreased, 
it evolves into the infrared (IR) region according to a couple of differential equations 



5.2 RG flow in the neighborhood of fixed points 

In the neighborhood of the respective fixed point {Xl,X2) in the plane {Xi,X2) = 
((^, A) or {a, a'), we can study the behavior of the RG fiow analytically. By taking 
into account only the terms which are linear in the infinitesimal deviation 6Xi := 
Xi — X^,6X2 := X2 — X2 from the fixed point, a set of RG equations, ( |5.5a| ) and 



( p.5b| ), is reduced to the form: ( j ~ A f j , where A is a two by two 



matrix. 





Eigenvalue 


Eigenvector 






(13 g^C2iG) . \ 


(e,A) 


{a, a') 




A 


1 


IR fixed point 


(l,a) 


any lines 


B 


-1 


UV fixed point 


(l,a) 


(l,a) 


any lines 


C 


1 


Saddle 


(13,3) 


(1,-2) 


line IV 




-1 


point 


(0,1) 


(0,1) 


line E 


D 


1 


Saddle 


(-13,3) 


(0,1) 


line V 




-1 


point 


(1,-2) 


(0,1) 


line n 



Table 1: Eigenvalues and eigenvectors of the linearized RG equation where the lines, 
I, n, IV are defined below. At the IR fixed point A and UV fixed point B, two 
eigenvalues are degenerate. 
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In A) plane, the set of linearized RG equations reads 



Similarly in (a, a') plane, we obtain 



3 J ' \la' I 3 (4 




The respective matrix characterizing the behavior of the RG flow in the neighborhood 
of the respective fixed point has the eigenvalue and the corresponding eigenvector 
enumerated in Table The direction of the flow is determined at the respective fixed 
point. We will see that these results are consistent with the global flow diagram given 
in Fig.^ below. 



5.3 Global behavior of the RG flow 

We find that ^ = 0,^ = 1/2, and ^ = 1 are solutions of the equation ( |5.5a| ). This 
implies that the RG flow starting from the point on one of the three planes, (0, A, g), 
{1/2, \,g), {l,\,g), is always kept on the respective plane. On the three planes, 
moreover, the equation ( p.5b| ) can be solved exactly. On the plane (1/2, A, g), the RG 
flow in the region < A < ^ obeys 




13 
22 



A(/x) = - 1+C In-^ , (5.11) 



where C is a positive constant. We see that A approaches to the ultraviolet (UV) 
fixed point A t y in the UV limit ^ oo, while A | monotonically as yU J. Aqcd- On 
the other hand, the RG flow in the region A > ^ is described by 



13 

22 



A(/i) = ? l-C^(ln-r^) } , (5.12) 



3 \ A 



QCD^ 
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where A approaches to the UV fixed point A t ^ in the UV hmit /i f oo, while A t C)0 



monotonically as | Aqcd- By substituting In ■ 



^QCD 



r 22 C2{G) 2\ 
\ 3 (47r)2 y / 



-1 



into the above 



equation, the equation of the RG fiow on the plane (1/2, A,^') is obtained 



l±C 



'22 C2{G) 2^ 
. 3 (47r)2 ^ ^ 



13 
22 



(5.13) 



The RG fiows on the plane (0, A, g) and (1, A, g) are governed by the same equations 
which are obtained by replacing 26/3 with 13/3. 

The global behavior of the RG fiow is obtained by solving ( 5.5c )-( ^.5b|) numeri- 
cally. In Fig.^, the RG fiow is drawn on the plane (^, A) and the plane (a, a'). The 
direction of the arrow denotes the direction towards the IR region and the length of 
the arrow is proportional to the magnitude of the vector /i^(^, A)/^'^. In the neigh- 
borhood of the respective fixed point, we see that the numerical result agrees with 
the analytical result given in Table |I] of the previous subsection. 



1 / 



n 



I': : 



I 



1 1 

4 » 
/ ( 
( 



\ \ 1 



i i 4 k k k k t / 
^ i i h k i k i / 



(b) 



A a' 



n 



< IB' 



, - - > I 



V 



a 



Figure 4: RG fiows in A) plane (a) and in (a, a') plane (b). 



Among the RG fiows, the five RG fiows (I, E, IE, IV, V) connecting the fixed points 
A,B,C,D form the watershed (or backbone) in the fiow diagram. 

a' = 0. (5.14a) 

a = 0. (5.14b) 

a' = -^a. (5.14c) 

a' = --a + —. (5.14d) 
2 3 

26 

a = -. (5.14e) 



I. 


e = 


2' 


E. 


e = 


0, 


EI. 


e = 


1, 


IV. 


A = 


13 1 
3 1 - 


V. 


A = 


131 

Ye 
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Since the flow is symmetric for the reflection with respect to the straight hne 
I:^ = 1/2, we focus on the region C, < 1/2. The flow starting from the initial position 
below IV runs towards the line A of fixed points and eventually arrive at A. If it arrive 
at a fixed point on A with a certain value of ^ depending on the initial position, then 
it does not move anymore. On the other hand, the flow starting from the initial 
position above IV runs away into the infinity, A = +00. Here the flow on the line I 
and E is not an exception. However, it should be remarked that the fixed point B 
is IR repulsive in both directions, while the fixed point C is IR attractive on IV and 
repulsive on I. In view of these, it turns out that any fixed point on A is IR stable, 
while the fixed point B on I is a rather special fixed point which is IR unstable (UV 
stable). g 

We have shown that the three fixed points B, C, D for the gauge parameter ^, A 
are located on the line I, E, III = 1/2,0, 1), respectively. On the lines I, E, H, the 
RG flow is confined in the respective line, the Lagrangian takes the following form. 

\. = \ (i.e., a G R, a' = 0): The GF+FP term is invariant under the FP ghost 
conjugation and the orthosymplectic transformation 05'p(4|2) [|T^] . 

£gf+fp = i8b5b {^A^ ■ - ■ e) . (5. 15) 

There is a four-ghost interaction. 

E. ,^ = (i.e., a = 0, a' G R): The GF+FP term is invariant under the global 
shift of anti-ghost C: 

a' 

^GF+FP = — S ■ S + S ■ a^^^ + iQ ■ d^'D^'[A]Q. (5.16) 

There is no 4-ghost interaction. This Lagrangian is the same as that in the 
conventional Lorentz gauge. 

n. ^ = 1 (i.e., a' = — 1«): The GF+FP term is invariant under the global shift of 
ghost C: 

ilQp^pp = ^3.3 + 3. d^A" + iQ ■ 'Df'[A]df,e. (5.17) 
There is no four-ghost interaction. 

The choice E or IE eliminates the four ghost interaction and the Yang-Mills theory 
reduces to the FP case. Once ^ = or ,^ = 1 is chosen, ^ is not renormalized by 
quantum corrections, since ^ = and = 1 are fixed point of the renormalization 
group. Then the FP Lagrangian is preserved under the renormalization. 

In E and EI, the role of ghost and anti-ghost is interchanged. The FP ghost 
conjugation invariance is broken in the usual FP Lagrangian where the ghost and 
anti-ghost are not treated on equal footing (except for the Landau gauge). In other 

^ This does not imply that the similar result is obtained also for the MA gauge. For example, 
a = is not a fixed point in the MA gauge. See j27) for details. 
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words, the FP ghost conjugation invariance is recovered for a' = (i.e., = 1/2 or 
A = a) by including the quartic ghost interaction even for a = 0. 

We must keep in mind that these results are obtained to one-loop order. Therefore, 
the details of the flow diagram may change if we include higher-order corrections. The 
higher-order result is not known to date and will be given elsewhere. Nevertheless, the 
existence of the fixed point at A = remains true to any finite order of perturbation. 
The existence of the lines, I, E and 111 are also guaranteed even after the inclusion of 
higher order terms, since it is protected by the symmetry dictated in the above. This 
is because the symmetry can not be broken as far as the perturbation series to all 
orders are not summed up. 



6 Renormalizing the composite operator of mass 
dimension 2 

In this section we discuss the renormalization of the composite operator of mass 
dimension 2 and its BRST and anti-BRST invariance under the renormalization. 



6.1 On-shell BRST transformation 



By eliminating the Nakanishi-Lautrup field the on-shell BRST and anti-BRST 
transformations are obtained as 

^-df^A^ix) + ^ige{x) x e(x) 



SBGix) 



A 



A 



d^Af.ix) - - l)ige{x) x e(x) 



(6.1) 
(6.2) 



The nilpotency of the on-shell transformation is partially broken]^ by the equation 
of motion of ghost and anti-ghost: 

{dsfA^ix) = (6.3a) 
0, (6.3b) 



(Ssmx) 



YM 



A 6e 

d^V^e - g^d^A^ X e) + tg^X^i^ - l)(e x C) x e. 



and 



(dBfA^ix) = 0, 
(^B)'e(a;) 

(^B)'e(x) 



xr tot 

T 5e 

d^D.e - g^d^A, X e) - ig'm - i)(e x C) x e, 

0. 



(6.3c) 
(6.4a) 



(6.4b) 
(6.4c) 



^ An elegant proof of the unitarity of the gauge theory is given based on the nilpotency of 
the BRST transformation, see e.g. ||3^. The nilpotency is indeed broken in the on-shell BRST 
transformation which is obtained by eliminating the NL field. However, the nilpotency is not the 
only way to show the unitarity. Even in this case, it is possible to show the unitarity order by order 
of the perturbation theory based on the Feynman diagrams without the NL fields. 
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Moreover, the anti-commutatibity is also broken in the similar way: 

{SbSb + SbSb)Ai,{x) = 0, 



{SbSb + SBSB)e{x) = 
{SbSb + SB6B)e{x) = 



1 Af t°t 

A 6e ' 

i OAjyM 



A 5e 



(6.5a) 
(6.5b) 

(6.5c) 



6.2 Composite operator of mass dimension 2 



Wc define the composite operator O as a linear combination of two composite oper- 
ators of mass dimension 2: 



O = (Q(^))-i J d^x -A^{x) ■ A^'ix) + >dt{x) ■ Q{x) 



The on-shell BRST transformation of the operator O is calculated as 



(6.6) 



-A^{x) ■ A^'ix) + \iQ{x) ■ Q{x) 



A^,[x) ■ 6bA^{x) - \it{x) ■ bB^{x) + \%8b^{x) ■ ^{x) 

A^{x) ■ d^eix) + Xie{x) ■ |(e(x) x e(x)) 



SbO ^ {n^''^)-^ J d^'x 5 
= (q(^))-i|a 

= (o(^))-iyci^x 

+ d^'A^ix) ■ e{x) - X^ig{e{x) x e(x)) ■ e{x 

= (Q(^))-^ J d^x [d^[A,ix) ■ eix)] + A - e) ■ 9ie{x) x e(x))} . 

(6.7) 

In the similar way, the on-shell anti-BRST transformation of the operator O is cal- 
culated as 

SbO = (n(^))-^ I d'^x |a''[yi^(x) • t{x)] + a - e) ^e(a;) • g{Q{x) X e(x))| . 

(6.8) 

Therefore, the composite operator O is invariant under the BRST and anti-BRST 
transformations when 



(6.9) 



i.e., on the line 1 and A in the A) plane, or on the line 1 in the (a, a') plane. For 
^ = 1/2, the on-shell BRST and anti-BRST transformation reads 



6bQ{x) = -'-d^'A.ix) - ]^gQ{x) x e(x), 

5bQ{.x) = +-d>'A^{x) - ]-gQ{x) x t{x). 
a I 



(6.10) 
(6.11) 
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The special case, A = (and a = to have a finite ^) is nothing but the Landau 
gauge in the conventional Lorentz gauge and the BRST and anti-BRST invariant 
operator O reduces to a simple form: 



->A^(a;) 



(6.12) 



Note that O' is invariant under the gauge transformation as well as the BRST and 
anti-BRST transformation. 



6.3 Renormalization of the composite operator 

Hereafter, we use the following notation to simplify the expressions. 

A:=A^, C:=e^, C := C^, 5 := S^. 



(6.13) 



We consider the Green function of the fundamental fields with an insertion of the 
composite operator of mass dimension 2. In the following, it is assumed that wc have 
already finished the renormalization for the fundamental field in the pcrturbativc the- 
ory. Therefore, we have only to consider the extra renormalization for the divergence 
coming from the inserted composite operators in the renormalized Green function. In 
order to take into account the operator mixing among composite operators with the 
same mass dimension and the same quantum number, we must introduce the matrix 
of renormalization factors Zx-, - ■ ■ , ^4: 



(6.14) 





'\AA 




( 








iCC 






Then, to the lowest nontrivial order, we find 



AA 

AA 

iCC 
iCC 



\AA 



+ 




+ 



iCC 



\AA 



iCC 



+ ^[ + y + 



+ V ) + ... 



where we have used the Feynman rule: 



(6.15a) 

(6.15b) 
(6.15c) 
(6.15d) 



(6.16a) 
(6.16b) 



with the dot denoting the insertion of a composite operator. 
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We show that the divergences coming from the compositeness are absorbed by 
taking the four renormahzation constants Zi, Z2, Z3, Z4 appropriately. The first 
example is 



AA 



\AA 



R 



Z^lAA 



\AA 



+ Z0IAA 



iCC 
+ 




+ 



Hence the lowest value of is 1: 



Zi = l + Z['^ + ■■■ . 
Zi (^iCC [\AA]) + Z2 {iCC[iCC] 



The second example is 



\AA 



R 



+ z, 



=0 (no divergence). 
Hence Z^ does not have the tree part and begins with the one-loop order: 



Z2 — Z^} + 



The third example is 



iCC 



iCC 



R 



Zs{iCC [lAA]) + Z4(iCC 



iCC 



z^ 

+ ^4 



Hence, Z^ has the form: 



The fourth example is 



Z4 = 1 + 4'^ + 



AA 



iCC 



R 



■■Z:,{AA 
Z^ 



\AA 



+ Z4 (AA 



iCC 




+ 



+ ^4 I V + + 

eO (no divergence). 



(6.17) 

(6.18) 



(6.19) 
(6.20) 



(6.21) 
(6.22) 



(6.23) 
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Hence, Z3 begins with the one-loop order: 

■^3 = -^3 ^ + 



(6.24) 



Therefore, up to one-loop order, the renormalization constants must satisfy the rela- 
tionship: 



(1) 




+ 



0, 
0, 



-r ^ r -r ^ y u. 



+ 



0. 



(6.25a) 
(6.25b) 
(6.25c) 
(6.25d) 



The explicit calculations lead to the following divergent part. 




3+-A(l + A) 



(47r)2 e' 



3 + A^(^^l 
4 (47r)2 e' 

e^2 1 



zC2(G')5^^e(l-0A'^^^ 



1 

-( 
4 



zC2(G)5^^e(l-OA 



(47r)2 e' 

(ff^l 
(47r)2 e' 

(47r)2 e' 
(^/x-)2 1 



(47r) 



(6.26) 
(6.27) 
(6.28) 
(6.29) 
(6.30) 
(6.31) 



Thus the renormalization constants for the composite operators are obtained as 



(1) 



--(1 + A)C.(G)-^^^^, 



= -A^e(l - 0C2(G) 



-c\2 



(47r)2 e' 



(6.32a) 

(6.32b) 

(6.32c) 
(6.32d) 



We pay attention to the renormalization constants of composite operators in light 
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of the inverted relation of ( |6.14| ): 



\AA 

\iCC] 



^3 



Z2 
^4 





'\AA 




)"'( 




:) 




[iCC 





(1) 




\AA 

\iCC] 



R 

R 



(6.33) 

This relation shows that there is an operator mixing between the gluon and ghost 
composite operators which are of mass dimension 2 and color singlet, as pointed out 
in [0]. In the absence of four-ghost interaction (^ = or ^ = 1), ( |6.28| ), ( |6.3CI| ) and 
( |6.31| ) vanish and hence we have = — '^^^'^ 



Z\ . In this case, there is no contribution 



from ghost for the renormalization of the gluon composite operator 

l'iCC]R - 



'-AA 

2 



'-AA 

2 



[iCC] 



7W 



R 

Ua 

2 



R 



(6.34) 
(6.35) 



On the other hand, the ghost composite operator can not be finite without the mixing 
of the gluon composite operator. In the conventional Lorentz gauge fixing, therefore, 
we do not have to consider the contribution from ghost in treating the renormalization 
of the gluon composite operator \AA (at least in the one- loop level). 



6.4 Multiplicative renormalizability of the composite opera- 
tor 

Now we examine the multiplicative renormalizability of the composite operator O. 
Taking into account the renormalization of the fundamental field and the composite 
field (|6.33|), we obtain 



Qo := 2 ^0^0 + AoiCoCo 



CC]^} 




-zx 

l + Z, 



\AA 



\{-Z, 



(1) 



\AA 



R 



z^-y^z!^) 



\AA 



+ 1 + Z« + Z« 



R 



zf\\\iCC\^,. 



(6.36) 



The multiplicative renormalizability holds (in the one-loop level) if and only if 



'(1) 



Z\ 



(1) 



\z. 



(1) 



zj_ 

A 



7(1) 



'(1) 



This is equivalent to the condition: 



0. 



(6.37) 



(6.38) 



26 



If this condition is satisfied, the composite operator is multiphcatively renormahzed 

as 



Qo 


= Zq{ 


7(1) 


/35 
~ V12 



^ + X[tCC]n), (6.39) 



-e\2 



In the case of A = 0, this result reduces to that of Boucaud et al. |^ without operator 
mixing. 

It should be remarked that the composite operator is not multiphcatively renor- 
malizable, unless the renormalization of the composite operators AA and CC are 
taken into account. In fact, the multiplicative renormalizability of 

Qo := ^AoAo + MCqCo 

1 + 4^) ^AA + (1 + 4) + 4')) XtCC + Oih''), (6.41) 

without the renormalization of the composite operator leads to the condition: 4^'' ^ 



Z^^ — 4^"* = 0, which reads A 



e(e 



|. This curve does not have a definite 



meaning in the renormalization, since the curve is not along the RG flow. 

6.5 BRST invariance of the renormalized composite operator 

Finally, we show that the renormalized composite operator is invariant under the 
renormalized BRST and anti-BRST transformation. By requiring that the renormal- 
ized BRST and anti-BRST transformations are nilpotent and anti-commute: 

(5M^0, 5|5| = 0, « + ^^(5^ = 0, (6.42) 

the renormalized BRST and anti-BRST transformation for the renormalized fields 
Afj^, C, C, B are determined (by an appropriate rescahng of 'B field) as 

5lA^{x) = XD^[A]^C(x) := X[d^C{x) + Z]l^Z,g^{A^{x) x C(x))], (6.43a) 
SlC{x) = -]^XZ)l^Zgg^{C{x) X C{x)), (6.43b) 

5lC{x) = tXB{x), (6.43c) 

d^B{x) = 0, (6.43d) 

and 

SiA^ix) = X'D^[A]''€ix) := X[d^Cix) + Z'f Z^g^iA^ix) x C'(x))], (6.44a) 

^B^(^) = -Ixz'fZggniCix) X Cix)l (6.44b) 

S^C{x) =iXB{x), (6.44c) 
S^B{x) = 0, (6.44d) 
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where X and X are arbitrary real numbers and B is defined by 

B{x) = ~B{x) + tZ]/^Zgg^{C{x) x C{x)). (6.45) 

The Lagrangian is written by making use of the renormalized BRST and anti- 
BRST transformation and the renormahzed fields as 

This agrees with ([4.21|) . 

We derive the condition for the renormalized composite operator O-^ to be in- 
variant under the renormalized BRST transformation defined above. We can write a 
finite composite operator of mass dimension 2 in the form (up to an overall constant): 



Q 



R 



\a,{x) ■ A^{x) 



+ K^[iC{x)-C{x)]n, (6.47) 

R 



where K^i is a finite function of the renormalized parameters, g-^, ^r, A. Performing the 
renormalized BRST transformation ( |6.43d| ) after the renormalization factors ( |6.33| ) 
of the composite operator are included, we obtain 

2r = Sl[{Zi + i^R^3)(^^^ ■ A^) + (Z2 + K^Z,){iC ■ C)| 
= (Zi + K^Z^)Xd^C ■ 

+ {Z2 + K^Z^)i^iC ■ {XZ^zf-C X C) 

+ X ^-^}^d,A^ - iZTz^iZ.gC X ■ C}. (6.48) 

For the right-hand-side to be a total derivative, we must require two conditions: 1) 
the coefficient for the term C ■ (C* x C) vanishes, 2) the remaining terms containing 
the derivative are combined into a total derivative term. The respective condition 
reads 



Z^ZgZ^i, (6.49) 

+ i^RZa = (Z2 + K^Z,)^^\. (6.50) 
The first condition reduces to 

io = Z^i = \. (6.51) 

Since Z2, ~ Oifije) and Zi, Z4 ~ 1 + 0{h/e), the second condition yields for the 
0(1) term: 



i^R = Ar, (6.52) 
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and for the 0(l/e) term: 

Zi^) - Z« - Ar4^) + ^ - - Zg) + Zl^) = 0. (6.53) 

This condition is the same as (|6.37|) . In the Landau gauge a = A = 0, especially, the 
condition (|6.53|) reduces to = 0. This is automatically satisfied in this case. 



7 Operator product expansion and vacuum con- 
densate 

We apply the operator product expansion (OPE) or short distance expansion (SDE) 
to the gluon and ghost propagators. The OPE is originally proposed as an operator 



relation by Wilson [56]. For example, the product of two scalar field operators defined 
at different spacetime points is expanded as 



(a;)0(y)~5]F['^'](x-i/) 



X + y 



(7.1) 



where the composite operators {Oi} form a complete set of renormalized local oper- 
ators. The famous proof of OPE by Zimmermann |^ was given in the framework 
of the perturbation theory. Quite recently, the OPE was rigorously proved as an 
operator relation by Bostelman [^j.Q According to the method p9|, pO|] , the (Fourier 
transformed) Wilson coefficient F['?^i '<^"](j9) in the OPE: 



(l){x)cl>{y)r~.J2F^^'-'^-\x-y) 



X + y 



(7.2) 



can be calculated in perturbation theory by equating a (2 + n)-point one-particle 
irreducible (IPI) Green's function — where two of the external legs have hard mo- 
mentum p and the remaining n external legs are assigned zero momentum q = — 
with the Wilson coefficient times an n-point Green's function with an insertion of the 
relevant composite operator at zero momentum. 



7.1 The OPE in the tree level 

First, we consider the OPE of the inverse gluon propagator: 

(D-yjip) = c^'^t^{p){i) + c^^'^t^{p){\A, . An + c^^'^^i^mc ■c) + ..., 

(7.3) 

where the first Wilson coefficient is nothing but the bare inverse gluon propagator: 



AB 



P^Y ^ ^-lP^\ ^AB 



p^ 



p^ 



(7.4) 



^ The authors would Ukc to thank Izumi Ojima for informing this reference. 
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The other Wilson coefficients are calculated in the perturbation theory from the 
diagrams: 




= ^ + ^ + O . (7.5) 




Ct'' = ^2\. + ■ (T.6) 



In the diagram, two external legs have hard momentum p and the {n = 2) lines 
connected to a blob correspond to the external legs with zero momentum g = 0. 

The explicit calculation in the tree level yields the result (see Appendix for the 
details of calculations): 

C^^^'i'ip) = 2(]|T^e(l - OPiJ^', (7.8) 

where we have put C2{G) = for simplicity. Defining the vacuum polarization 
tensor of the gluon by 

(D-X^'^ip) := (D.X'^ip) + X^^ip), (7.9) 
we obtain the vacuum polarization tensor of the gluon: 

K^ip) = 4(^2!^ 1) ^^'^ {-(1 + ^)P^^(^p ■ ^0 + 2^^(l - OPiui^C -C)}. (7.10) 

It turns out that even the inclusion of the quartic ghost interaction does not affect 
the Wilson coefficient Cj^^^ at least in the tree level. For the Wilson coefficient Cj^*"' , 
however, there is an extra contribution coming from the quartic ghost interaction, as 
suggested already in 0]. The non-zero Wilson coefficient C^'"'"^ due to the presence of 
the quartic ghost interaction (^ 7^ 0, 1) breaks the transversality of the gluon polariza- 
tion tensor, i.e., U^^ 7^ Pju^- This result does not contradict with the Slavnov- Taylor 
identity When ^ = (resp. ^ = 1), the ghost condensate {iC ■ C) can not 



appear in the OPE, since the gluon-ghost-anti-ghost vertex ( [4.16| ) is proportional 
to the outgoing ghost (resp. anti-ghost) momentum p^ (resp. q^). The above result 
( [7.1U| ) suggests the existence of the effective gluon mass given by 

< = -^^^^(l + ^)(^p-^0- (7.11) 

Therefore, the gluon condensation of mass dimension 2 can be an origin of the gluon 
mass. The effect of higher orders will be investigated in the next subsection. 
Next, we perform the OPE for the inverse ghost propagator: 

-^{G-Y''{P) = 4'kp)(l) + C^aI\p){\a, ■ A") + cf^\pWC ■ C) + ■ • ■ , (7.12) 
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where the first Wilson coefficient agrees with the bare inverse ghost propagator: 



= -KG^y'ip) = (7.13) 
The other Wilson coefficients are calculated from the diagrams, 

^ -a, 



which yield the result: 



C]S?(P) = 0. (7.17) 

\CC] 

Here the coefficient vanishes due to cancellation, see Appendix. Defining the 

vacuum polarization tensor of the ghost by 

(G-Y^'ip) ■■= (Go ^)^^(P) + ^Ili^{p), (7.18) 
the vacuum polarization for the ghost is obtained: 

K'(p) = ^^^f'i^p-^')- (7.19) 

We find that the ghost vacuum polarization has no contribution from the ghost-anti- 
ghost condensation even for ^ 7^ 0, 1. Thus we obtain the effective ghost mass: 

This result shows that the gluon condensation of mass dimension 2 can also be an 
origin of the ghost mass.^ 

The combination of gluon and ghost condensation appearing in the OPE is not 
BRST invariant in the sense explained in the previous section. This is reasonable, 
since even the OPE of gauge invariant operators does not give a gauge invariant 



combination in the OPE, see e.g. |51 



®In the Lorenz gauge, the effective gluon mass and ghost mass are generated by the gluon conden- 
sation of mass dimension 2 alone in the tree level. This is not the case if we include the high-order 
correction as will be shown in the next subsection. In the MA gauge, on the contrary, two conden- 
sations from the off-diagonal gluon and off-diagonal ghost contribute to the effective off-diagonal 
gluon and ghost masses already in the tree level, see P,[2^. 
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7.2 RG improvement of the OPE 

One of the advantages of the OPE is that the momentum dependence of the Wilson co- 
efficient is determined by the renormahzation group (RG) equation. More accurately, 
the change of the Wilson coefficient under the RG transformation can be specified 
by the renormahzation factors Z which are to be calculated before the RG improve- 
ment of the OPE calculus. Therefore, we can obtain the higher-order corrections 
for the momentum dependence of the coefficient without any explicit higher-order 
computations (at least for the leading logarithmic corrections). 

7.2.1 RG equation for Wilson coefficients 

We begin with an OPE relation in the momentum representation obtained by ex- 
tracting composite operators up to mass dimension 2 (we omit all the indices, since 
they are not essential in the following arguments): 



iA(0)A(0) 



+ F2^(p)[zC(0)C(0)]r + -- - , 

R 

(7.21a) 



1a(0)A(0) 



^ + F,^{p)[tC{0)Cmn + --- , 

(7.21b) 



where Dp^,t{p) and Gpert denote respectively the perturbative gluon and ghost prop- 
agators in which the perturbative loop corrections are included besides the OPE 
contribution. 

First, we try to rewrite all the field operators in both sides of (|7.21a|) and (|7.21b|) 



in terms of bare quantities. Hereafter it is supposed that the Wilson coefficient and 
composite operators are defined based on the renormahzation scheme depending on 
a certain parameter fi (corresponding to the mass scale), which is different from 
the BPHZ prescription at zero momentum, g = 0. In the actual calculations, we 
adopt the minimal subtraction (MS) scheme, although the resulting expressions can 
be translated into those of the momentum-space subtraction scheme (MOM). 

By making use of the Z factors calculated already in the previous section, two 
OPE relations above are combined into a matrix form: 

1 / -zio(p)io(-p) \ _ / iAo(0)Ao(0) \ 

' [ doip)Coi-p) ) - + [ ^Co{0)Co{0) j + " " " ' ^''^^^ 

where we have defined the two by two matrices, 

y _( Za \ ^._( F^{p) Fiip) \ ~_(Z,Z,\ , 



and a column vector, 
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Introducing a matrix Fq by | 

Fq = ZfFZ : = 



Fotip) Fo^ip) 



(7.25) 



we obtain an OPE relation among the bare quantities as 



Oo{p)Coi-p) 



7n A-V ( Mo(0)^o(0) 



+ 



(7.26) 



Second, we observe that the relation (|7.26 ) should have no dependence on the 
renormalization point /i. Hence, the first term on the right- hand- side of ( [7 .261 ) is 
independent of /x, i.e.. 



/x^(ZfDp,,J = 0, 

and the coefficient Fq in the second term is also independent of /i, i.e, 

/i— Fo =/i— (ZfFZ) = 0. 
a/i ct/i 



(7.27) 



(7.2J 



We multiply ( [7.28|) by Z^ ^ from the left and by Z ^ from the right to obtain 



where denotes the parameters of the theory (5'r,^r, Ar), and denotes the cor- 
responding RG function Pi{a) := /i^ttj- Here we have used a fact that /i^ -|- 
Y.i is just the ordinary differential operator //^. 



Defining the RG function (matrix) 7f,7 from Zf, Z by 

d ^ rv d ~ _ ~ 

/i— Zf := Zf7f, /i— Z := 7Z, 
d/i djjL 



(7.30) 



we obtain the RG equation for the matrix F of the Wilson coefficients: 

F(p, a, /i) + 7fF(p, a, /i) + F(p, a, /i)7 = 
Similarly, we can show that Dp^rt obeys the RG equation: 

Dp,rt(p, a, /i) 7fDp,rt(p, a, /i) = 0. 



(7.31) 



(7.32) 



''Were it not for the renormalization of the composite operator, Fo reduced to F. 
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7.2.2 Solving the RG equation 



Now we proceed to solve the RG equation just obtained. A simple dimensional 
analysis leads to a relation, F{K,p,a, nfi) = k'^^F {p,a, fi) which is equivalent to the 
relation: 



F{kp, a, /i) = k'^^F [p,a,— 



/i 



(7.33) 



d d 

+ 

on oil 



(7.34) 



where (if is the canonical dimension of F. Hence, F satisfies 

F(kp, a, /i) = 0. 
We use this equation to eliminate jj,-^ in ( [7.31|) to obtain 

F(kp, a, yu) — 7fF(Kp, a, /x) — F(/tp, a, /i)7 = . (7.35) 



da. 



This is the homogeneous RG equation of Weinberg-'t Hooft type which is adequate 
for the mass-independent renormalization method. 

By the standard method ||3y,^, the general solution of the RG equation ( [7.35| ) 
is given by 



F {K,p,a, fj,) = K,~^ exp <. J dn' — - — > F(p, «(«), /x) exp < y dK 



(7.36) 



where we have imposed the boundary condition: a{K, = 1) = a. 

The similar consideration yields the general solution of the RG equation ( [7.321 ): 



Dp.rt {Kp,a,fi) = K ^exp 



dn' 



Dpcrt (p, a, /u) • 



(7.37) 



Once we know the Z factors of the fundamental field and of the composite op- 
erator, it is easy to calculate 7f, 7 according to ( [7.30|) . If the integrations in the 
arguments of the exponential in ( |7.36|) and (|7.37| ) are performed, the k dependence of 
the solution will be exactly determined. However, Z factors are obtained in terms of 
renormalized parameters g-^, ^r, Ar and hence depend implicitly on k through them. 
This fact makes the analysis more difficult in general. 



7.2.3 Solution around the UV fixed point B 

We can calculate 7f, 7 up to 0{h), since we have known all the Z factors of the 
fundamental field and of the composite operator up to 0{h) in this paper. In the 
high-energy limit k — 00, it is expected that the solution can be explicitly obtained 
in the neighborhood of the non-trivial UV stable fixed point in the parameter space, 
due to asymptotic freedom of the Yang-Mills theory, i.e., q^k) — = as k ^ 00. 

In the three-dimensional parameter space ^r, Ar, actually, we have found that 
all the points are fiowing into the UV fixed point B in the UV limit except for some 
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lines that have higher symmetry. On the other hand, within the perturbation theory 
using the dimensional regularization, the /i dependent loop correction of all the Z 
factors always appears with a factor of Oi^g^. Therefore, the RG function 7 as an 
element of the matrix 7 defined by differentiating the Z factor with respect to [i 
is accompanied by g\ to the 0{K), like 7 ~ S'r/I^, + Oif?^. If the polynomial 
function /(^, A) in the above expression has a nonvanishing value at the fixed point 
(^*,A*) , the \i dependence of 7 = g^j is governed by alone and hence we can 
replace /(^, A) with the constant /(^*, A*) at the UV fixed point. By substituting the 
fixed-point values A^ = 26/3, = 1/2 into ^, A, the Z factors become 

13^?2iV,/i-2^ YJg'N^y.-^^ 
^ 4 167r2 e ' ^ 4 V 3 / 167r2 e ' 

which yield the matrix of the renormalization group function: 



gj, T-W^l^^ |. (7.39) 

Furthermore, we define the coefficient matrix C^j and in (|7.39|) by 

7;(^?):=/C,„ ^\g):=g'C^. (7.40) 

By taking into account the RG equation ^ij^g = —-^g^ (b = ^^c) and the result- 
ing relation: ■^Ing'^ = ^-^g = the nontrivial integration of { 7.3t ) can be 
performed as 



Hence the solution becomes 

F{kp, a, /i) = K-^ { ^ 1 (F(p, «(«:), /i)) ( ^ ) . (7.42) 
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.g'^{K) ' ' V5'^('«J 



The K dependence of g"^ is obtained by solving its RG equation as ^^(k) ~ In /t]^^ 
for large k. Substituting ( [7.41j ) into (|7.36|) , therefore, we determine the In/t depen- 
dence of the solution for large k: 

F(Kp,a,/i) = /t-^(lnK)*^^f^ (F(p,a(K),/i)) (In/t)*^^^ . (7.43) 

In order to cast the matrix power of In/t into a more tractable form, we shall 
diagonalize the matrix C^y in such a way that S diagonalizes by the similarity 
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transformation —>■ S ^C^S. Such a matrix S and the diagonahzed matrix are 
given by 

/_13 _26\ N n\ 

I). (7.44, 

This diagonahzation corresponds to redefining the combination between two com- 
posite operators of mass dimension 2, i.e., ^A{0)A{0) and iC{0)C{0), by multiplying 

3_( IA^ + fiCC \ 
= 13 - f^Cc) ■ 





(7.45) 

Inserting the identity matrix 1 = SS^^ appropriately, the solution ( |7.42|) is rewrit- 



ten as 





F{kp, a, /i) = K-"" ^ F(p, a, ij)SS~' ^ SS''. (7.46) 



Now both C^^ and S ^C^yS are diagonal. Hence we can write down the power explic- 
itly as 



13 Nc \ / o 3 iV, 



F(kp) = K 



^4 



U^WJ T(p)S ^^^W^ 

I (SS)^^^v V 



4 2b 








(7.47) 



Here we impose a condition that T(p) := F(p, «(/€), /i) coincides with the Wilson 
coefficient in the tree level obtained in the previous section in which the coupling 
constant is replaced with the running coupling constant a{K). Note that F is the 
Wilson coefficient of the Green function (not of the one-particle irreducible (IPI) 
function) .[^ Hence we put 



T(p) 



(Trip) T2{p)\ 
\Up) Up)) 

N,-g\K) (-{iD,)\l + X)P^ (zDo)'4e(l-Oa 



2(iV2 - 1) \^ {iG^)' 



(7.48) 



We notice that each element Ti,-- - ,T4 of T(p) brings an extra Iuk factor to F 
through (7^(k) ~ Therefore, the OPE correction up to dimension 2 operators 
reads 

F(p) 




3 -r J2nin^/AQCD'' 3 ^1 ^ ^2nin^/AQCD'' | 

137^ /Inp/AQCDNM:^ 267^ / Inp/Agcp \^-^ \ 



(7.49) 



Except for the Landau gauge in which any operator mixing does not occur, a linear combina- 
tion of different powers of In k appears in the solution, and its combination coefficients cannot be 
completely determined by perturbation theory alone. But it is important to note that a fitting of 
the analytical result with the simulation data (or experimental data) can determine the asymptotic 
behavior of F completely as discussed in the next subsection. 
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where we have used T4 — 0. Here we have used the translation rule from the MS 
scheme to the MOM scheme: 



Inp/AqcD 



g'^{K) In/z/AqcD' 



(7.50) 



Among the terms with various powers of In k, the largest-power term (correspond- 
ing to the largest eigenvalue of the matrix C^) is dominant in the UV limit (k ^ 1). 
Extracting this In k, contribution, we can simplify the Wilson coefficient of IPI func- 
tion in the UV limit as 



C 



ipi 



gh '-'gh 



-2 



pert / 





2b (N^ - 1] 










(Gpert/ Go 



-2 



\-2 



(7.51) 



X 



13(-l-A)PT-6$(l-g)PL ' AqCD 



(In ■ 



13 



^ ^QCD 



(In- 

26 Nc 



'^QCD 



61 Nn , 

61 'Vc 
)T2-W 



13(-l-A)PT-6^(l-g)PL -^QCD 



(In- 



'QCD 



26 Nc 

)"B""5r 



13 



^ ^QCD 



(In- 
26 Nc 



^QCD 



61 Nc , \ 

61 A'<: 
)T2 2b 



QCD 



26 Nc 



(7.52) 



In the similar way, we obtain 



Dpert(«P,a;,/^) = « 








17 iVc 
12 26 



Dt(p), 



(7.53) 



where the tree expression is given by 
Dt(p) - 



Doip) 
-^Go(p), 



-^(Pt + APl) 



(7.54) 



7.2.4 The solution at the conventional Landau gauge 

Finally, we consider the OPE on the line A of the fixed points (corresponding to the 
conventional Landau gauge), the RG matrices read 



7f 



7f 



87r2 I 







7 



87r2 



12 
_i 

2 



0^ 

4y 



(7.55) 



The diagonalization can be performed as 



'0 -f- 



87r2 lo i, 
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(7.56) 



The eigenvalues of C-^, are the same as those at the fixed point B. Therefore, we 
obtain the Wilson coefficient C]£'^ between {A^{p)A^{-p))-'^ and ((^(0))^) and C^^^"^ 
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between {C{p)C{-p))~^ and ((^(O))^): 



3 iVe 

4 2b 



Y{kp) 



K 



where no mixing between gluon and ghost occurs due to T2 
The coefficient of the IPI OPE read 



(7.57) 

in addition to T4 = 0. 



a 



2b 2{N^ - 1) 



D 



pert 



,(ln. 



Do 



'QCO 



3 JVc 



pert ^ 



,(ln- 



^ '^QCD ' 

13 Nc -I 



Go 



'QCD 



where 



Dpert(p) 



Mnp/AqcD y 





13 Nc 
6 2b 



(In 







'^QCD 



13 Wc 



(7.5^ 



/ Inp/AQCp V 
Vln/x/AqcD/ 



iv, D,(p). 



3 iV^ 

4 2b 



(7.59) 



This result for the ghost part is new, while the gluon part reproduces the recent 
result of Boucaud et. al. in the MOM scheme (Note that their definition of 7 is 
different from ours by a factor 2 and the coefficient 70 differs by the signature). In 
order to transfer from our renormalization scheme to the MOM scheme, we have used 
the translation rule ( [7.5U| ). In the Landau gauge, therefore, we have conffimed that 
the ghost condensation does not affect the inverse gluon propagator as in the tree 
level, even if the leading logarithmic corrections are taken into account in the OPE. 
In other words, the gluon condensation is decoupled from the ghost condensation 
within this approximation. 



7.3 Full propagators: momentum dependence 

The vacuum polarization tensor of the gluon is decomposed into the transverse and 
the longitudinal parts: 

K'^ip) = [n^(j>')^^. + n^(i>')^^J5^^ (7.60) 

where n""" and II'" are functions of alone. Once the vacuum polarization functions 
n""" and n'" of the gluon are obtained from the OPE, the propagator is written as 



-p2 + nT(p2) 



^ 

X 



where we have defined a function Z^i{—p^) by 

Zg\{-P^) = Zp,,t(-p^) + Zope(-P^) := 



XU^ip 



P.. 



,2^ ^ Ml' 



-P 



-p2 + nT(p2)- 



(7.61) 

(7.62) 
(7.63) 
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Note that 11^ (p^) = in the conventional Landau gauge. 

On the other hand, if the vacuum polarization function of the ghost Ilgjf (p^) 
S'^^Ugti^p^) is calculated by OPE, the ghost propagator is obtained as 

ip -\- liigiiyp ) p 

where we have introduced a function Ggt{—p'^) by 



(7.64) 



V) = G,U-P') + GoM-p') ■■= _ ( (7-65) 

P + '-^g'ayP 



The OPE contribution Yl^^^ to the vacuum polarization function in the inverse 
propagators ( |773| ) and ( [7. 121 ) is related to the Wilson coefficient C^^^ as 

Substituting the result ( [7.49| ) into (|7.66|) , we obtain a pair of vacuum polarization 
functions: 



'(rj^ 1371 y lnp/AQCD Nf^:^ 1 ( rj. 26 m W In p/Agcp N 61 :^ 1 \ . 



(7.67) 

It turns out that the vacuum polarization functions just obtained reduce to the tree 
results, i.e., ( |7.1CI| ) and ( [7.19| ), at k = 1 (or p = /i). Therefore, the ghost condensation 
(iCC) does contribute to the gluon and ghost vacuum polarization functions in the 
leading logarithmic corrections of the OPE. 

Thus the OPE contribution to the gluon and ghost vacuum polarization functions 
are obtained: 

35 Nc -I 

_2vr2iV,(l + A)/(lnA^)--- / 1 26 



12 2b 



i (In-^) 



^ V2 ' ■ 3 



61_iVc_-i , 



(7.68) 
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(7.69) 
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The effective gluon mass is obtained from the pole of Zgi(— p^), i.e., a solution of 
the equation = UtIp^), while the effective ghost mass is obtained from the pole 
of Ggh(— p^), i.e., a solution of the equation = —illgh{p'^). In view of this, the 
solutions ( |7.68| ) and ( [7.69| ) would give an improvement of the tree- level result, ( [7.11|) 



and ( [7.20| ). However, a BRST non-invariant combination Q2 of composite operators 
appears together with the BRST invariant combination Qi discussed in the previous 
section. Therefore, these results indicate that we need more endeavor in order to 
reach the BRST invariant pole position in the IR region. 
In the Landau gauge, especially, we have 



3 iV, 




4 2b ^ 



ZA-P') = ~P'DUP)-P-' I Y (N2 u 3^ W I , (7.70) 

13 Nc -f 

Ggh(V) = -^p'GUp)+p-' I Y r/V2 n MiY. {A') 1 . (7.71) 

After the Wick rotation to the Euclidean region p"^ — > —p%, we find that the function 
Zgi(p|;) is monotonically increasing in p|; if := — (v4^) > 0, as in the case of 

constant U^{-pl) = > 0. On the other hand, if := -{A^) < 0, ^gi(p|) 

has a Landau pole in the IR region and is monotonically decreasing in p"^ in the 
UV region. In the conventional Landau gauge, these results can be compared with 
those of the Schwinger-Dyson equation (see e.g., |5^) and the numerical simulation 
on a lattice (see e.g., P^|5^ -|5B|). According to these results, Z^p^) is enhanced at 



intermediate momenta and has a peak at about 1 GeV. It was argued [^] that the 



enhancement of the gluonic form factor at IR region is related to quark confinement. 
However, this region is beyond the reach of our study in this paper. Incidentally, the 
data in the gauge other than the Landau gauge is not yet available. 



8 Conclusion and discussion 

In this paper we have discussed the multiplicative renormalizability of the composite 
operator O in QED and Yang-Mills theory. This research is motivated by clarifying 
the mechanism of mass generation and a possible connection to quark confinement. 

In QED, we have shown that the composite operator is trivially renormalizable 
and that the renormalized composite operator is BRST and anti-BRST invariant for 
an arbitrary value of the gauge fixing parameter. There is no subtlety related to the 
renormalization of the composite operator. 

In the Yang-Mills theory, we have adopted the most general Lorentz gauge with 
two gauge- fixing parameters ^, A which was derived by the Baulieu and Thierry- 
Mieg IQ]. It was known that the bare composite operator O of mass dimension 
2 is invariant under the hare BRST and anti-BRST transformation for the choice 
of gauge parameters A = or ,^ = ^ and that it is also invariant under the gauge 
transformation in the Landau gauge A = 0. In this paper the composite operator has 
been renormalized by taking into account the operator mixing carefully. Here the 
matrix of renormalization factors has been explicitly calculated. Consequently, we 
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have found that the BRST and anti-BRST invariance of the renormahzed composite 
operator holds if the renormahzed parameters take the same value, Ar = or 
■Cr = |) the bare one. Moreover, we have obtained the KG flow in the (^, A) plane 
to one-loop order. In the RG flow diagram, the RG flow runs only on the line = | 
if the initial position of ^ is located somewhere on the line. The line Ar = is a 
line of fixed points. Therefore, if the system is located on a point in the line Ar = 
initially, it can not move from the initial position. This fact guarantees the BRST 
invariance of the renormahzed composite operator O^. 

We have also examined in this paper how the conventional calculations are mod- 
ified in the presence of the vacuum condensate of mass dimension 2. By performing 
the OPE of the gluon and ghost propagators, we have shown that the effective masses 
of gluon and ghost are generated due to the non-vanishing vacuum condensate. Al- 
though this phenomenon was already suggested based on the tree level calculation, 
we have taken into account the leading logarithmic corrections in consistent with the 
RG flow by making use of the RG equation. We have found that the effective masses 
are provided from the ghost condensation {iQ ■ C) as well as the gluon condensation 
{^A^ ■ A^) (except for the Landau gauge A = 0). This result should be compared 
with the tree level result where the effective mass has the contribution from the gluon 
condensate alone. 

The next step is to show that the non-vanishing vacuum condensates (O) ^ 
is actually realized in the QCD vacuum. An attempt in this direction has already 
been performed in |^ by calculating the effective potential for the ghost conden- 
sation (iCC) in the SU{2) and 5*^7(3) Yang-Mills theories in the MA gauge. Quite 
recently, Verschelde et al. have carefully obtained the multiplicatively renormal- 
izable effective potential for the gluon condensate Q^A^A^) in the Landau gauge up 
to two-loop order in the SU{N) Yang-Mills theory. Both results support that the 
non-zero vacuum condensate of mass dimension 2 is energetically favoured in the 
Yang-Mills theory. In these approaches, an auxiliary field p{x) corresponding to the 
composite operator has been introduced to obtain the effective potential V{a) of a 
constant a = p{x). However, this treatment has a number of subtle points which 
have not been discussed in these papers. This issue will be discussed in a subsequent 
paper |]2^ in detail. 

In massless QED, photon pairing |^,^ can occur in the strong coupling phase 
||39|-^ where the chiral symmetry is spontaneously broken. Therefore, it will be 
possible to discuss the interplay between quark confinement and chiral symmetry 
breaking on equal footing in a unified treatment. The extension of this viewpoint 
into the non-Abelian case, i.e., gluon pairing |^ is also an interesting subject to be 
tackled in the future work. 

Finally, we point out that the operator O is essentially a mass term for the gluon 
and the ghost fields. Although a naive introduction of a mass term for the gluon 
alone breaks the BRST symmetry, our result indicates that there is a BRST invariant 
combination of mass terms: 

"1 

This mass term is very similar to that obtained after the spontaneous breakdown 
caused by the non-vanishing vacuum expectation value of the Higgs scalar field. In 



■ — tr 



-m^Af,{x) ■ A^{x) + m^aie{x) ■ e{x) 
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our case, the mass should be of dynamical origin. It is possible to give a proof 
of the multiplicative renormalizability of the Yang-Mills theory with a mass term 
preserving the BRST symmetry to all orders of perturbation theory. However, it is 



known [^, ^ that the introduction of the mass term breaks the nilpotency 

of the off-shell BRST transformation as well as the on-shell one. Consequently, the 
unitarity of the theory turns out to be spoiled. In this sense, the mass generation 
should occur in the dynamical way, i.e., {O) 7^ in the limit m — * 0. This viewpoint 
will be discussed in a subsequent paper. 
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A OPE calculations 

In order to give the OPE correction for the gluon propagator, we need to calculate 
the following diagrams. 



^ :w\-N;)^[{A + X)p'Pl + {D- l)p'PlJ\5^^. (A.l) 




(iV- 

+ f^"" f '''''' {9,.9,8 - 9,89,.) + f^'' f '''''' {9,u9,8 - 9,,98u) 
■■'9 \^,l^^j^ N^[9,u{D - l)]5^^. (A.2) 



2 



(A.3) 



(iV2-l)p2' 

For the correction of the ghost propagator, we need the calculation of the following 
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diagrams. 



2 



(iV2 - 1) 
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